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Abstract 

A method for solving Schwinger-Dyson equations for the Green 
function generating functional of non-Abelian gauge theory is pro- 
posed. The method is based on an approximation of Schwinger-Dyson 
equations by exactly soluble equations. For the SU(2) model the first 
step equations of the iteration scheme are solved which define a gauge 
field propagator. Apart from the usual perturbative solution, a non- 
perturbative solution is found, which corresponds to the spontaneous 
symmetry breaking and obeys infrared finite behaviour of the propa- 
gator. 



Introduction 



The Schwinger-Dyson equations (SDE) method is one of the basic 
tools for investigating of the Green functions of the quantum the- 
ory. Hitherto the one universal method for solving SDEs has been the 
coupling constant perturbation theory (hereafter simply termed per- 
turbation theory). The range of applicability of other methods (for 
example, 1/iV-expansion) is limited by a narrow class of models. In 
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particular, the 1 /^-expansion method cannot be used in the calcula- 
tions for non-Abelian gauge theories due to the complicated structure 
of the leading approximation. 

On the other hand, the applicability of perturbation theory to 
the investigation of non-Abelian gauge theories is limited by a deep- 
Euclidean region. In the non-perturbative region of small momenta 
the physical vacuum of non-Abelian gauge theories obeys the nontriv- 
ial structure that is beyond the framework of perturbation theory. In 
the SDE terms this fact can be understood if one takes into account 
the radical difference in the properties of the leading-approximation 
equations of the perturbation theory and the original exact equations. 
The SDEs for the generating functional of Green functions are equa- 
tions in functional derivatives. The leading approximation of per- 
turbation theory comprises neglecting the higher derivatives terms in 
these equations (just such terms correspond to an interaction). The 
leading-approximation equations of perturbation theory are of a lower 
order compared with the exact ones; therefore, the class of solutions 
thus described contracts drastically, and non-perturbative solutions 
which correspond to the nontrivial physical vacuum practically drop 
out of the consideration. This feature of SDEs in the non-perturbative 
region is noted repeatedly in simple models (see Q). 

In this work a method for the SDE solution of non-Abelian gauge 
theory is proposed. It takes into account the terms with higher deriva- 
tives (i.e. self-interaction of the non-Abelian fields) ab ovo in the 
leading approximation. Though we limit ourselves to the simplest ex- 
tension of the class of SDE solutions, the results are non-trivial: the 
non-perturbative solution which corresponds to spontaneous symme- 
try breaking (with non-Higgs mechanism) is found. It obeys a non- 
singular behaviour of the propagator in the non-perturbative infrared 
region of small momenta. 

An idea of the method comprises the approximation of SDEs for 
the generating functional by equations with "constant" (i.e. indepen- 
dent of sources) coefficients. These equations have a simple exponen- 
tial solution, which is a foundation for the linear iteration scheme. 
The method is universal, as is perturbation theory, i.e. it is applica- 
ble to practically any model of quantum field theory. For the scalar 
(f) A theory it has been shown in that the method describes such 
non-perturbative phenomena as spontaneous symmetry breaking and 
the trivialization of c/> 4 theory at d = 4. This method has also been 
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successfully applied in the investigation of the Gross-Neveu model for 
finite N §]. 



2 Schwinger-Dyson equations and the 
iteration scheme 

A system of SDEs for the generating functional G{J,rj) of the Green 
functions of non-Abelian gauge theory has the form 



^<(7T7)^ + ^ = 0. (2) 



Hdj'Sff 

Here F£ V (A) = d^A a u - d v A° + gf ahc A\A% is the gauge field tensor, 
T>^(A) = S^dp—gf^A*^ is the covariant derivative, f abc are structure 
constants of a gauge group, J^(x) is the source of the gauge field, r] a (x) 
is the source of a ghost field, a is the gauge parameter, g is the coupling 
constant. We work in Minkowski space with a metric (1, -1, -1, -1), 
and = g^ u x^y u by definition. 

At g = the system (l)-(2) has a solution 

Gf ree = exp{^r J * Df ree * J + ifj * A free * 7]}, 

where Df ree and Af ree are free propagators of the gauge field and 
ghost field respectively, and 

J * D fr ee * J = J dxdyJ^x)(D free ) a ^(x - y)J^y), 

etc. The functional Gf ree is the generating functional of the free Green 
functions. The iteration scheme based on Gf ree is the usual coupling 
constant perturbation theory. 

We shall use an alternative iteration scheme for SDEs (l)-(2), 
which is formulated as follows: the leading approximation is a system 
of equations with terms proportional to the sources J and rj omitted. 
(For the system of eqs. (1) and (2) these are the last terms.) This 
system has a solution 

Go = expi{J *V + fj*C + C*7]}, 
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where J * V = f d x J®(x) V®(x), etc. The coefficient functions 
Vj% and C a are solutions of the corresponding system of " characteristic 
equations" . When constructing the iteration scheme for the generating 
functional 



the omitted terms JG and rjG should be considered as perturbations, 
i.e. equations of the iteration scheme are 



A solution of eqs.(3)-(4) has the form G n = P n Go, where P n is a 
polynomial in J and rj. Therefore, at each iteration step we obtain 
a closed system of equations for coefficient functions of the polyno- 
mial P n , which completely defines the Green functions of the given 
step. There is no manifest small parameter in the usual sense in this 
scheme: "smallness" is defined by the condition that the Green func- 
tions are derivatives of the generating functional at J = r\ = 0, and it 
is sufficient for us to know G{ J, rj) near zero, i.e. in the region where 
the neglected terms are small. At each step of the iteration we ap- 
proximate the functional G/Gq by a sum of the polynomials P n , and 
a degree of the polynomial increases with each step. As is known, for 
ordinary differential equations the scheme of this type is equivalent to 
iterations of Volterra-type integral equations and gives an expansion 
that converge well. That is why we may hope that this scheme pos- 
sesses good convergence properties. In any case it is clear that the 
convergence of this scheme is no worse than that of the perturbation 
theory. The perturbation theory is singular in the sense of differential 
equation theory, since the higher derivatives are omitted in the lead- 
ing approximation. In contrast to the singular perturbation theory, 
the scheme proposed is regular in the sense above. This circumstance 
gives us hope for improving the convergence properties (see [§] for 
more discussion). 

To remove the ultraviolet divergences it is necessary to supple- 
ment SDEs (l)-(2) and the iteration scheme equations (3)-(4) with 
the corresponding counterterms. The counterterms are also defined 
by the iteration procedure: Sz = Szq + 5z\ + • • •, i.e. at each step it is 



G — Gq + G\ + • • • + G n + • • • 



J/j,G n —i, (3) 
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necessary to take into account the counterterms of the corresponding 
order. 

Let us consider the leading approximation in more detail. As has 
been noted above the solution of the leading approximation equations 
is the linear exponential in the sources. The vacuum structure is 
defined by the solutions of the characteristic equations for and C a 

V V (V)F VI1 (V) + -dpdvVy + gfCd^C = 0, (5) 
a 

d^(v)c = o. 

These equations have a great number of solutions which reflects the 
non-trivial vacuum structure of the non-Abelian gauge theory. In this 
connection it is appropriate to recall 't Hooft's conjecture pi on the 
existence of different vacuum modes for non-Abelian gauge theory 
at zero temperature: a superconducting one (spontaneously broken), 
dual superconducting (confinement), etc. Realization of either mode 
depends on the values of some quantities (manager parameters) whose 
definition is a problem of the quantum field dynamics. It is reasonable 
to suppose that the whole set of solutions of characteristic equations 
(5) is divided into classes, each corresponds to some mode of the gauge 
theory in the above sense . 

In this paper we limit ourselves to investigation of the first step of 
the iteration scheme based on the simplest solutions of the character- 
istic equations (5). 

The trivial solution of eqs.(5) V£ = C a = (i.e., G = 1) leads 
to an iteration scheme which is simply the reconstructed perturba- 
tion theory. The simplest extension of the set of vacuum solutions 
compared with the trivial perturbative vacuum comprises taking into 
account non-zero constant solutions of eqs.(5). Namely, we choose the 
leading approximation in the form 

G = expi{J*V}, (6) 

where the vacuum vectors V£ are independent of the space-time vari- 
ables and satisfy the condition 

fab c f cdhybydyh = ^ (?) 

which follows from the characteristic equations (5), if C = dV = 0. 
Below we shall consider only this class of the leading approximation 
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solutions. Surely we have no prior physical foundation (except for 
reasons of simplicity) for the choice of this class of solutions. Never- 
theless, we shall see that this simple set of solutions leads to a new 
spontaneously broken mode. 

Concluding our discussion of the leading approximation, let us 
touch upon the definition of the ground state ("physical vacuum"). 
Due to the linearity of SDEs (l)-(2) and the iteration scheme equations 
(3)- (4) an arbitrary linear combination of solutions is also a solution. 
Therefore, the physical vacuum functional G =< | >j should be 
constructed as a linear combination of partial solutions, each of them 
corresponding to some solution of the characteristic equations (5). 
In other words, the physical vacuum is a superposition of F-vacua. 
This combination should be chosen in such a way as to maintain the 
admissible physical properties of the Green functions (for example, 
< | A a | >= 0) and the energy minimality condition. 

To formulate the iteration scheme in terms of the polynomials P n , 
it is convenient to introduce the matrix quantity 

W ab = igf abcyc (g) 

Then, the iteration scheme equations have the form of equations for 
the polynomials P n : 

m{ v) + */j*r][iW JL) + w.±. - w,±] + 

w ^-k + +g, k4 }p " = - J " p "- 1 ' (9) 

dMV)+isf-£-]? P f=i'lPn-l- (10) 

Here Pq = 1. The first-step solution Pi defines the leading approxima- 
tion for the propagators of the gauge and ghost fields. The following 
steps define the many-particle functions. 



3 Solution of first-step equations 

The solution of first-step equations of the iteration scheme (9)- (10) is 
the quadratic polynomial in the sources 

P 1 = —J * D * J + if] * A*r]. (11) 
2z 
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Eqs. (9) and (10) give us equations for the functions D^ v {x — y) and 
A ab (x — y). The equation for can be essentially simplified by 
modifying the gauge condition. Instead of the usual covariant gauge 
used above, it is convenient to use the modified gauge condition ("V- 
gauge" ) with the gauge fixing term 

C aaug e = ~{V ll {V)A li ) 2 . (12) 

The ghost terms should be changed correspondingly. For the transi- 
tion to the V-gs,xige in formulae (l)-(4), (5) and (9)-(10), it is sufficient 
to perform the substitution 

d^V^V). (13) 

An essential circumstance is that the leading approximation condition 
(7) is not changed in the case. 

The equations have a particularly simple form in the gauge a = 1 
("diagonal F-gauge"). Then, the equation for D^ u in momentum 
space is 

{/CV + 2[K li ,K v ]}D vX {k) = -<7 MA , (14) 

where k is the momentum, and the following notation has been intro- 
duced: 

/qf = k^ ah - Wf. (15) 

In the region of large k eq. (14) tends to the equation for the free 
propagator in the diagonal gauge, i.e. at k — ► oo 

Dpvik) « ~~tf.9p.v- ( 16 ) 

The equation for the ghost propagator A is 

K, 2 A(k) = -1. (17) 

In the large k region the propagator A also tends to the free prop- 
agator. Therefore, the ultraviolet behaviour of the solutions with a 
nontrivial vacuum vector V M is the same as for the usual perturbation 
theory. 

Below we consider the case of SU(2) gauge group and restrict 
ourselves to the subset of constant with zero field tensor F^ U (V) = 
0, or, equivalently, with 

e abcybyc = Q _ ( lg ) 
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It is clear that condition (18) ensures the leading approximation con- 
dition (7). The energy of that kind of solutions is zero, as for the 
perturbative solution with trivial vacuum. In that sense this set of so- 
lutions can be termed quasiperturbative. Eq. (18) gives us [!C^,IC U } = 
[Wp, W u ] = 0, and the solution of eq.(14) is reduced to the inversion 
of the matrix K? . The result is 

^ pi pi k pi V z 

where the following notations has been introduced: 

p{k, V) = k 2 + g 2 V 2 , Pl {k, V) = p 2 - 4g 2 {kV) 2 . (20) 

A solution of the ghost propagator equation (17) is also given by 
formula (19) (without g^)- 

Except for eqs. (14) and (17), the first-step equations give one 
more relation that contains a quantity D^ u (0), which should be un- 
derstood as some regularization. In essence this relation is a condition 
for the first step counterterms bz\. (There is no need to introduce the 
leading approximation counterterms in this case, i.e. 5zq = 0). Since 
the Green functions of the first step are finite, this condition for the 
counterterms is necessary to remove the ultraviolet divergences at the 
second step of the scheme. This peculiarity of the given iteration 
scheme is displayed here in exactly the same manner as for the scalar 
field theory (see Q). 



4 Gauge field propagator 

Let us turn now to a possible physical interpretation of the solutions. 
At V M = D^ u and A are free propagators of gauge and ghost field, 
and the whole iteration scheme is a reconstructed series of the pertur- 
bation theory. At V M ^ the situation is more complicated. It is clear 
that in this case it is difficult to interpret the function given by 
(19) as the propagator of a particle in the Poincare invariant theory. 
Let us recall, however, that we have a number {G(V)} of solutions of 
SDEs, each corresponding to some vector satisfying the condition 
(18). Due to eq.(18) the field strength tensor is identically zero and 
any such solution possesses zero energy. Therefore, a candidate for 
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the "physical vacuum functional" for this set of solutions is a super- 
position of all solutions. We shall exploit this fact in the construction 
of a Poincare invariant solution which can be interpreted as a particle 
propagator, i.e. a function depending only on the momentum k and 
the scalar quantity 

v 2 = v 2 = v;v;. (21) 

(The quantity v 2 plays the role of an order parameter.) The construc- 
tion is equivalent, in essence, to some averaging, i.e. an integration 
with a measure dfj,(V), we shall therefore denote it by angular brack- 
ets: < G >= J2vG(V), etc. In the foundation of this operation we 
set the following conditions 

< y« >= 0, < V 2 >= v 2 . (22) 

Their necessity for the Poincare invariant theory is evident. It is also 
evident that 

< V*V u b >= \v 2 g, u E a \ (23) 

where tr E = 1. To determine the form of the matrix E ab , consider 
the leading approximation condition (18). The geometrical meaning 
of condition (18) is the collinearity of the vectors V M in isotopic space. 
Consequently, at ^ there exists a selected direction in isotopic 
space. This direction can be chosen as a basis vector, for example n 3 . 
In this basis V£ = <5 a3 ^, and E ab = 5 a3 5 b3 . Therefore, at V M / the 
isotopic symmetry is spontaneously broken. 

Further calculation is reduced to the averaging of the functions 
f({kV) 2 ) = p/pi and f ab = f ■ (V a V b )/V 2 . First of all, note that at 
k — > / — > l/g 2 v 2 , and at k — > oo / — > l/k 2 . These Poincare- 
invariant properties of / should, of course, be conserved for < / > as 
well. 

With eqs.(22)-(24) and its generalizations for an arbitrary mono- 
mial in Vu, the following formulae can be proved 

< <^ 2 " >= ' tV '~ ^p - < 24 > 

< (W) 2n {V a V b ) >= E ab v 2 < {kV) 2n >, 

which are necessary for the calculation of < D^ v > . The result of the 
calculation in the above basis is as follows: 

< Bf u {k) >= -g^, (25) 
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/f)H^W /f>22,,w g + ffV A /! 4fl 2 t> 2 fc 2 \ 

(26) 

(Other isotopic components are equal to zero.) Therefore, along the 
selected isotopic direction the particle propagates as a free one, but 
along other directions a separation of the region of momenta exists, 
the scale of the separation being the quantity g 2 v 2 . It is necessary 
to stress that both the limiting cases above (k 2 — > and k 2 — > oo) 
belong to the region of applicability of the calculations performed, 
which is defined by the condition | Ag 2 v 2 k 2 / (k 2 + g 2 v 2 ) 2 |< 1. Con- 
sider the question about analytical continuation. Eq. (26) defines two 
analytical functions depending on the choice of branch of the func- 
tion V^i 2 , but neither satisfies both the above asymptotical conditions 
simultaneously and, consequently, these functions are not solutions 
of the problem. Hence, as a solution one should choose at large k 
the branch with the behaviour 1/A; 2 , and at small k - another branch 
which is a constant l/g 2 v 2 . At the points ±g 2 v 2 the solution goes 
from one branch to another, i.e. from the " perturbative" sheet to 
the "non-perturbative" one. At the point k 2 = g 2 v 2 the solution is 
continuous, and at the point k 2 = —g 2 v 2 has a discontinuity. (Note 
that the first-step calculations of the iteration scheme do not fix the 
sign of v 2 .) Such unusual features of the solution near the separation 
points are likely to be connected with our limitation of the characteris- 
tic equation solutions. Probably an extension of the class of solutions 
will lead to smoothing of the propagator behaviour near the points of 
separation of the perturbative and non-perturbative regions. 



5 Discussion 

In this paper we have attempted to introduce a new iteration scheme 
to solve the Schwinger-Dyson equations for non-Abelian gauge the- 
ory. As a very first step we have considered a scheme based on the 
simplest constant solutions of the characteristic equations (5). Of 
course, there are no prior grounds for expecting such a set of solu- 
tions to be useful for the physical applications. Some other sets of 
more non-trivial solutions (instanton like and so on) of the charac- 
teristic equations are likely to be of real physical interest. Dealing 
with coordinate-dependent solutions is a more difficult task; in partic- 
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ular, the averaging procedure of Sec. 4 should be essentially modified. 
Another important question is the correspondence of the proposed 
method with the functional integral approaches, such as well-known 
background field method ||. In spite of their evident similarity, the 
proposed iteration scheme seems to be somewhat different: the struc- 
ture of expansions is distinct. The problem of the correspondence 
requires further investigations based on the calculation of higher steps 
of the iteration scheme. 

The set of quasiperturbative solutions of Sects. 3-4 leads to spon- 
taneous breaking of SU(2) gauge symmetry. The question arises as 
to whether such solutions can be used to construct a realistic model 
of electroweak interactions in the spirit of dynamical gauge symme- 
try breaking models without Higgs bosons (see @ and refs. therein). 
Further investigation which should essentially extend the class of so- 
lutions is necessary to eliminate this question. On the other hand, the 
set of quasiperturbative solutions satisfying condition (18) considered 
here does not exhaust the solutions with zero energy, and some further 
averaging to restore the symmetry may be necessary. 

For SU(3) gauge group (the case of QCD) the variety of solutions 
increases greatly in comparison with SU{2) group (see, for example, 
[||) . Avoiding discussion of this case, nevertheless note that the 
infrared finite behaviour of the gluon field propagator has recently 
been discussed in detail (see and refs. therein). 
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